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Abstract —In this paper, we consider a MU-MISO system where 
users have highly accurate Channel State Information (CSI), 
while the Base Station (BS) has partial CSI consisting of an 
imperfect channel estimate and statistical knowledge of the CSI 
error. With the objective of maximizing the Average Snm Rate 
(ASR) subject to a power constraint, a special transmission 
scheme is considered where the BS transmits a common symbol in 
a multicast fashion, in addition to the conventional private sym¬ 
bols. This scheme is termed Joint Multicasting and Broadcasting 
(JMB). The ASR problem is transformed into an angmented 
Average Weighted Sum Mean Square Error (AWSMSE) prob¬ 
lem which is solved using Alternating Optimization (AO). The 
enhanced rate performance accompanied with the Incorporation 
of the mniticast part is demonstrated through simulations. 

Index Terms —Joint Multicasting and Broadcasting (JMB), 
Imperfect CSIT, Robust Design, AWMSE. 

I. Introduction 

The availability of accurate Channel State Information at 
the Transmitter (CSIT) is crucial for Downlink (DL) Multi- 
User (MU) multi-antenna wireless transmission. This stems 
from the necessity to deal with the interference through 
preprocessing at the transmitter side, as the receivers are 
distributed 12 . While the ability to provide highly accurate and 
up-to-date CSIT remains questionable, considerable effort has 
been devoted to improving the performance in the presence of 
CSIT uncertainties. Recent information theoretic developments 
focusing on the Multiple Input Single Output (MISO) Broad¬ 
cast Channel (BC) suggest that multicast assisted transmission 
(where common symbols which are decodable by all users are 
transmitted alongside the conventional private symbols) can 
be used to enhance the performance in the infinite Signal to 
Noise Ratio (SNR) regime El, El- This paper focuses on the 
particular case where linear precoding is employed to transmit 
one common symbol in addition to private symbols in each 
channel use. The simultaneous utilization of the MU-MISO 
medium as a Multicast Channel (MC) lIU and a BC ||5] is 
termed Joint Multicasting and Broadcasting (JMB). 

For a CSIT error that decays with increased SNR, JMB 
was shown to boost the achievable sum Degrees of Freedom 
(DoF) m, 0. However, these results are intrinsically focused 
on the asymptotic SNR regime, where the DoF analysis is 
most meaningful. Therefore, trivial choices of linear precoders 
are deemed sufficient given that they achieve the aspired sum 
DoF. For example, naive Zero Forcing Beamforming (ZF-BF) 
is used for the private symbols, while the multicast precoder 


is not optimized. However, this is not the case at finite SNR 
where more involved performance metrics are considered, 
e.g. the Sum Rate (SR). Works that consider the BC and 
MC separately under simpler CSI assumptions (i.e. perfect 
CSI) suggest that the choice of precoders can significantly 
influence the performance E), G]. However, the instantaneous 
SR cannot be considered as a design metric at the BS due to 
the CSI uncertainty. Alternatively, the Averag^H SR (ASR) is 
considered as an overall performance metric. ASR maximiza¬ 
tion problems are tackled by extending the approach in jS], 
i.e. transforming them into augmented Average Weighted Sum 
Mean Square Error (AWSMSE) problems which are solved 
using Alternating Optimization (AO) 0, noi. 

Contribution and Organization: In this work, we employ 
JMB transmission to optimize the SR performance for a MU- 
MISO system with partial CSIT. To the best of our knowledge, 
this has not been considered in literature. Due to the stochastic 
nature of the CSIT uncertainty, precoders are designed such 
that the ASR is maximized. The problem is transformed into 
an equivalent augmented AWSMSE problem, solved using an 
AO algorithm which converges to a stationary point. Moreover, 
we demonstrate the benefits of incorporating the common 
symbol. In particular, it is shown that the asymptotic DoE gains 
translate into SR gains in the finitely high SNR regime. On 
the other hand, JMB reduces to conventional MU transmission 
when the common symbol is not needed, e.g. at low SNRs. 

The rest of the paper is organized as follows: the system 
model and problem formulation are described in Section HI] 
In Section HII] the equivalent augmented AWSMSE problem 
is introduced. An AO algorithm that solves the AWSMSE 
problem is proposed in Section |IV] Simulation results are 
given in Section |V] and Section |Vl| concludes the paper. 

Notation: Boldface uppercase letters denote matrices, bold¬ 
face lowercase letters denote column vectors and standard 
letters denote scalars. The superscrips (•)^ and (•)^ denote 
transpose and conjugate-transpose (Hermitian) operators, re¬ 
spectively. tr(-) and || • || are the trace and Euclidian norm 
operators, respectively. Einally, Ea;{-} denotes the expectation 
w.r.t the random variable x. 

II. System Model and Problem Eormulation 

We consider a Base Station (BS) equipped with Nt antennas 
serving K {K < Nt) single-antenna users. The BS operates in 

* The term ’’Average” is used to denote the expectation w.r.t the CSIT eiTor. 


a JMB fashion transmitting K private symbols, each intended 
solely for one user, in addition to a common symbol that is 
decodable by all users. The vector of complex data symbols 
is given as s = [sc, si,..., Sir]^ G where Sc is the 

common symbol, Si is the ith user private symbol, i G 1C and 
/C = {1,..., K}. Entries of s have zero-means, unity powers 
and are mutually uncorrelated such that E{ss^} = I. s is 
linearly precoded into the transmit vector x G given as 

K 

X = PcSc + ^PiSj (1) 

where pc G C^* and pi G C^* correspond to the precoders for 
the common symbol and the jth private symbol respectively, 
from which P = [pc, pi,..., p/f] is composed. The total 
transmit power at the BS is denoted as Pt, from which the 
transmit power constraint is given as E{x^x} = tr(PP^) < 
Pt- For the fcth user, the received signal is given as 

Uk = hf X + Uk (2) 

where G is the narrow-band channel impulse response 
vector between the BS and the fcth user, from which the 
composite channel is defined as H = [hi,...,hx]. Uk ~ 
CA/"(0,ct^j,) is the AWGN at the fcth receiver with variance 
. Throughout the paper, it is assumed that noise variances 
are equal across all users i.e. = <7^; \/k G IC. 

A. CSIT Uncertainty 

Each of the K links exhibits independent fading, and 
remains almost constant over a frame of symbols, enabling 
users to estimate their channel vectors with high accuracy. 
On the other hand, CSIT experiences uncertainty arising from 
limited feedback, delays or mismatches. H is written as a sum 
of the transmitter-side channel estimate H = [hi,..., h^f] 
and the channel estimation error H = [hi,..., h^f], such 
that H = H -f H. The CSIT consists of H, in addition to 
some statistical knowledge of H. Particularly, the BS knows 
the probability distribution of the actual channel given the 
available estimate, i.e. /h|h(H) = /h(H — H). 

B. MSB, MMSE and Rate 

The kth user obtains an estimate of the common symbol by 
applying a scalar equalizer pc,fc(hfc) to (|2| such that 'Sc,k = 
9c,k0^k)yk- Assuming that the common symbol is successfully 
decoded by all users, the common symbol’s receive signal 
part is reconstructed and cancelled from yk- This improves 
the detectability of Sk, which is then estimated by applying 
5 fc(hfe) such that Sk = 5 fe(hfc)(yfe - hf PcSc.fe). The notations 
5 c.fc(hfc) and gk{irk) are used to emphasise the dependencies 
on the actual channel, as each user is assumed to have perfect 
knowledge of its own channel vector, (h^) is omitted for 
brevity unless special emphasis is necessary. This is used with 
other variables that depend on the actual channel. For the fcth 
user, the MSEs defined as ec,k = Es,rafc{|sc,fc — Scp} and 
£k = Es,„fc{|sfe - Sfcl^} are given as 

ec.fe(hfe) = |5c.feP'Ec,/c - 2K{pc.fehf Pc} -f 1 (3a) 

£k(.hk) = IpfcpTfe - 2ifi{gkh^Pk} + 1 (3b) 


where T^^k = jpf -f Tk and Tk = |pf h/cP + ul- 

The optimum pc,/c and gk are obtained by setting and 

to zeros, yielding the well known MMSE equalizers: 

9c,k (n/cj-Pc and pfc (hfej-PfchfcJ^ . (4) 

Substituting (|4|i into Q, the fcth user’s MMSEs are given as 

e^^^^{hk) = T-^E,^k and ef^^^{hk) = T-^Ek (5) 

where Ec,k = Tc,k - jpf = Tk and Ek=Tk- jpf h^l^. 

The MMSE and the Signal to Interference plus Noise Ratio 
(SINR) are related such that yc,k = (1 — and 

= (1 _ gMMSEy^MMSE^ where yc,k and yk are the fcth 
user’s SINRs. Therefore, the fcth user’s maximum achievable 
common rate and private rate are written as i?c,fe(hfe) = 
— log 2 (e^^®^) and RkO^k) = — log 2 (£^'^®®), respectively. 
The common message is transmitted at a common rate defined 
as Rc = m\nj{Rc,j}f^i, which ensures that it is decodable 
by all users. Ultimately, the objective would be to design P 
such that the SR given as Rc + EfeLi E maximized. It 
can be seen that in scenarios where the multicast part is not 
beneficial, allocating zero power to the common precoder will 
yield R^ — 0, and the system reduces to a conventional BC. 
However, rates are functions of the actual channel and hence 
cannot be used to construct an optimization problem at the BS. 
Alternatively, we consider the Average Rates (ARs) defined as: 
^h,\h,{Rc,k} and Ei^^|£ji?fc}. In the following, Ej^^|gJ-} 
will be simply referred to as E{-}. Before we proceed to 
the ASR problem formulation, we highlight the benefit of 
incorporating the common symbol from a DoF perspective. 

C. DoF Motivated Design 

The DoF-motivated JMB design in El is briefly revisited 
in this subsection. Consider cr^ = 1 ^ SNR = Pt, and 
an average estimation error power E{||hfe||^} that decays as 
where a > 0 is an exponent that represents the 
CSIT quality. For example, a = 0 represents a fixed error 
power w.r.t SNR, e.g. constant number of feedback bits. On 
the other hand, a = oo corresponds to perfect CSIT. In DoF 
analysis, it is customary to truncate the exponent such that 
a < 1, where a = 1 corresponds to perfect CSIT from a 
DoF perspective |j2l. Under these assumptions, the precoders 
of the private symbols are given as = yjP^ /Kpjf, 
where p^ is a normalized ZF-BF vector constructed using the 
channel estimate H, such that ||pf^|| = 1 and h.fp^ = 0, 
\/i,k G IC, i ^ k. The common symbol’s precoder is given as 
Pj. = ^Pt — P^Bi, where ei is a standard unity basis vector 
with 1 as the first entry and zeros elsewhere. The fcth user’s 
received signal is given as 

yk = hf PcSc + hf pfcSfc -f ^ hf piSi -f 

o{Pl) o(pg) o(p«) O(p0) 

indicating the order of the average power of each term as 
Pt 00 . The third Right Hand Side (RHS) term corresponds 
to the residual interference from unintended private symbols. 





resulting from the employment of an imperfect channel esti¬ 
mate to construct the ZF-BF vectors. Since the error scales as 
O whilst the power allocated to the private precoders 

scales as O {P^), the residual interference is drowned by noise 
and can be neglected. By decoding the common symbol while 
treating the rest of the terms as noise, a DoF of 1 — a is 
achieved. Moreover, the private symbol achieves a DoF of a 
after cancelling the common symbol. The same applies to the 
other users, and a sum DoF of 1 + {K — 1 )q; is achieved. On 
the other hand, excluding the common symbol and splitting Pt 
between the private symbols, the receive power of the intended 
private symbol is enhanced to 0{P^). However, residual 
interference is also increased to 0(P/““), and the sum DoF 
obtained by the private symbols remains as Ka. Therefore, 
JMB is strictly superior to ZF-BF and SU transmission (e.g. 
TDMA which achieves a DoF of 1) for 0 < a < 1. The reader 
is referred to El, El for more on DoF analysis. It is clear 
that the DoF-motivated scheme adapts to the CSIT accuracy 
by changing the power allocation. However, the fact that the 
precoders are not optimized leaves considerable potential for 
improvement, particularly in the finite SNR regime. 


objective is to design P that maximizes the ASR defined as 
Rc + J2k=i subject to a power constraint Pt. This problem 
is formulated as 

K 


TZ- : max 
flc.P 

i?c + Rk 

(7a) 

k—1 


s.t. 

Rc,k > Rc, Vfc e /C 

(7b) 


tr(PP^) < Pt 

(7c) 


where the constraints in (TTbl ) are introduced to eliminate the 
potential non-smoothness arising from the pointwise mini¬ 
mization in Re- Problem 7?. is a non-convex optimization 
problem that appears to be very challenging to solve. 

HI. AWSMSE Optimization 
In this section, the ASR problem is transformed into an 
equivalent problem that can be solved using AO. We start 
by introducing the main components used to construct the 
equivalent problem, i.e. the augmented WMSEs jSl : 

Cc.fc(hfc,Hc.ic,tic.fc) — tic,fc(hfc)£^c,/c(hfc) log 2 (h/c)) (8a) 
Cfc (h/i;, (8b) 


D. Problem Formulation 


In order to formulate a deterministic ASR problem, the 
stochastic ARs are approximated by corresponding Sample 
Average Functions (SAFs). Each SAE is obtained by taking the 
ensemble average over a sample of M independent identically 
distributed (i.i.d) realizations drawn from the distribution /jj|h 
in a Monte-Carlo fashion. The sample is defined as Hm = 


{H(^ 


&M], where = [h 




, is the 


mth realization, and M = 


as: Rl^k 

where i?' 


M 2 -^m—l 


a... 


, M}. The SAEs are given 


and R 


(M) 


(m) 


— — U 

M Z^m=l 




■c.k = a,fc(hP) and R]^> ^ i?fc(hr ) are the 
rates associated with the realization . In the following, the 
superscript (m) is used to indicate the association of variables 
with the mth Monte-Carlo realization. It should be noted that 
P is fixed over the M realizations of the rates, which follows 
from the definition of the ARs. This also reflects the fact that 
P is optimized at the BS using partial CSI knowledge. 


Assumption 1. In the following, we assume that > 0 

with probability \,\lk € 1C. 


Alternatively, we can say that SNR = Ptlcr^ can only grow 
finitely large, and channel gains are finite. Assumption[T]yields 
> 0 with probability 1, as the presence of a 
nonzero noise variance dictates that Ec,k,Ek > 0. This also 
implies that rates are finite, and by the strong law of large 
numbers we can write 

Rc.k — lim = E{i?c.fc}, almost surely (6a) 

M—)-00 ’ 

Rk = lim = E{i?fc}, almost surely (6b) 

M—^CO 

where R^^k and Rk are the approximated ARs for a sufficiently 
large M, which will be simply referred to as the ARs. 
The common AR is defined as R^ = minj{Rc,j}jLi- The 


where Uc,fc(hfe) > 0 and Ufc(hfc) > 0 are weights associated 
with the fcth user’s MSEs, and the dependencies in (l8]l are 
highlighted for their significance in the following analysis. The 
dependencies of the weights on the actual channel is crucial for 
the establishment of the following WMSE-Rate relationship: 

min = 1 - Rc,k and min 4 = 1 - i?fc. (9) 

'^c,k •<9c,k ^kiQk 


This can be shown as follows: from = 0 and ^ = 0, 

ogc.fc agk 

the optimum equalizers are obtained as p* ^. = and 

g* = Substituting this back into ®, we obtain the 

augmented WMMSEs written as 

yilk^Uc^kj — ^c,k^c,k mg2(tlc,fc) 


Sc.fc 

tMMSE/v \ ^MMSE , „ / n 

{nk,Uk)=Ukek -\og 2 [Uk). 


(10a) 

(10b) 


Eurthermore, from 


q,MMSE 

duc.k 


= 0 and 


cMMSE 


the optimum MMSE weights: u*f. = u^ 


and ut = u 


MMSE 


= 0, we obtain 
A ^gMMSE^-1 

= Substituting this back into 

([Tol l yields the relationship in (|9]l. It is evident from (|5]l that 
the MMSE weights are dependent on the channel. 

The equivalent problem is formulated using the augmented 
AWMSEs defined as: E{^c,fc} and E{^fc}. Before we proceed, 
the augmented AWMSEs are approximated as: 


M 


FiM) ^ 

Z_^ Sc,, 


9 


M 


and 


c(M) _ J_ 

jvf Sfc 


9 


where 


& = & »S') and «<"' 4 4 . (h<”‘! „(”l) 


correspond to the mth realization of the augmented WM¬ 
SEs, which depend on the mth realization of the equalizers: 

dc^k - and = gk(h^^), and the weights: 

Uc^^ — and = Uk(h^^\ Eor compactness, 

we define the set of equalizers associated with the M real¬ 

izations and the K users as: G = {gc,fc,gfe | fc S A9}, where 






gc,fc = \ m G M} and gfc = { 5 ^™^ | m G Ad}. In a 

similar manner, we define: U = {uc^fc,Ufc | fc G /C}, where 

Uc,fc = I m G M} and | m G Ad}. The 

approximated augmented AWMSEs for a sufficiently large M 
are defined as Cc,fc = limM^oo^cX' Cfe - 
which will be simply referred to as the AWMSEs. The same 
approach used to prove ® can be employed to show that 
min Cc.fc = 1 - -Rc./c and min = 1 - .Rfc ( 11 ) 

^CjfciSc.fc UfciSfc 

where optimality conditions are checked separately for each 
realization. The sets of optimum MMSE equalizers associated 
with (fTTT i are defined as | m G Ad} 

and g^^®^ = I ^ g j\/[j jji the same man¬ 

ner, the sets of optimum MMSE weights are defined as 

MMSE(m) ^ r MMSE(m) , ^ A MMSE(m) ^ 

^c.fc — Wc,fc I ^ ^ — 

I^MMSE(m) I ^ g users, the MMSE solution 

is composed as | fc G /C} and 

UMMSE a {uMMSE^ „mmse I ^ ^ 

A. Augmented AWSMSE Minimization 

Motivated by the relationship in (HUl, the augmented 
AWSMSE minimization problem is formulated as 

K 


min 

Cc+ 

( 12 a) 

^■c,P,U,G 

fc=r 


s.t. 

ic,k < 'ik GK, 

( 12 b) 


tr(PP^) < Pt 

( 12 c) 


where is an auxiliary variable. The MMSE solution of the 
problems in (fTTl i is not only optimum for problem but it is 
also at the heart of a relationship that connects the stationary 
points of problem A to the stationary points of problem IZ. In 
the following, the notations and u'^'^®^(P) are 

used to emphasize the dependencies on particular precoders. 

Proposition 1. Eor any stationary point of A given as 

P*, U*, G*) that achieves an objective function of 
there exists: 1) a corresponding MMSE stationary point 
given as p*, u““®®(P*), G“^®J^(P*)) that achieves 
the same objective function, 2) a stationary point of TZ given 
as (1 — P*) that achieves an ASR of AT + 1 — Einally, 

if (f*, P*, U*, G*) is a global optimal point of A, then 
(1 — P*) must be a global optimal point of 'JZ. 

This can be shown by employing the ideas used to prove 
M Proposition 1]. A sketch of the proof is given as follows. 

Proof of Proposition\I}i Erom the KKT conditions of A^ 
it can be seen that (g^^®^, u^^®®) is optimal and unique 
Vk G iC. Moreover, (g^^®^, u^^®^) is optimal and unique 
Vk G /Ca, where /Ca C is the set of active constraints 
in (fT2bll . Eor inactive constraints, the MMSE solution is not 
unique but it satisfies the optimality conditions. This proves 
the first part. The second part is proved by examining the 
KKT conditions of problem 7Z and employing the relationship 
in The final part is proved by contradiction. Eor the 

complete proof, readers are referred to the extended version 
of this paper ca. ■ 


IV. Alternating Optimization Algorithm 
A lthough problem A is non-convex in the joint set of 
optimization variables, it is convex in each of the blocks P, U 
and G, assuming that the other two are fixed. This block-wise 
convexity is exploited using an AO algorithm that switches 
between optimizing blocks. Each iteration of the proposed 
algorithm consists of two steps: 1) updating G and U for 
a given P, 2) updating P for given G and U. 


A. Updating the Equalizers and Weights 

In nth iteration of the AO algorithm, the equalizers and 
weights are updated such that G = G’^“®®(P) and U = 
pjMMSE^p^ respectively, where P is the precoding matrix 
obtained in the (n — l)th iteration. To facilitate the problem 
formulation in the following step, the AWMSEs are written 
in terms of the updated blocks G and U, and the block P 
which is yet to be update. Eor this purpose, we introduce 
the AWMMSE-components listed as: ^k, tc,k, ik, tc,k, 

ffe, Uc,k, Uk, Vc,k and Vk, which are obtained using the 
updated G and U. In particular, the components Uc^k and 
Uk are calculated by taking the ensemble averages over the 


M realizations of and 


0 


. The rest of the components 
are calculated in a similar manner by averaging over their 
corresponding realizations given as: 


Am) 

^c,k 

(m) 

= 

(m) 

9 c,k 

2 

and 

II 

I.. 

■40 


= 4"A 

(m), 

i “fc 

and 

II 

Ig. 

Am) 

^c,k 

= «S’| 

(m) (m)^ 

‘■k 9c,k 

and 

II 

(m) 

^c,k 

^ logs 

[if 


and 

(m) 

^fe = 




(m) 

9 k 


,(m), (m), (m 
7c '■'■k '■'■k 

(m), (m) (m) 

^k “fe 9k 

")■ 


M 


H 


l 0 g 2 (t 


The AWMSEs are written in terms of the updated (G, U), 
and P which is yet to be updated, as 

K 

fc,fe = pf ^c.fePc + X! pf’^c.fcPc + cr^fc.fc - 2 K{fi}fcPc} 

H“ (133.) 

K 

Cfe ='^pf'^kPi +crltk - 2Ji{f^pfc} +Uk-Vk. (13b) 

i=r 


B. Updating the Precoders 

Eollowing the previous step, the problem of updating P 
is denoted by Ap, which is formulated by substituting (fTSI) 
into (fT^ and eliminating (G, U) from the set of optimization 
variables. This is given as 


K 


Cc.P 


K 


Ap :min Cc + X! ( '^Pf'^kPt - 25i{f^pfe} j (14a) 


k—1 ^i—1 
K 


s.t. pf’5'c.fcPc+X]p*^'®'cAP*+c^«^c.fe-23?{fi;^fePc} 


2=1 


A tXcjfc Pcjfc A ^k G KL 
tr( PP^) < Pt 


(14b) 

(14c) 


where the constant term + Uk — Vk) has been 

omitted from (fl4a . Problem (fl4l i is a convex Quadratically 






Constrained Quadratic Program (QCQP) which can be solved 
using interior-point methods ifTSll . 

C. Alternating Optimization Algorithm 

The AO algorithm is constructed by repeating the steps 
described in the two previous subsections until convergence. 
This is summarized in Algorithm [T] where en determines the 
accuracy of the solution and rimax is the maximum number of 
iterations. Step [T] is discussed in Section IVl 


Algorithm 1 Alternating Optimization 

1: Initialize: n ^ 0, ^ 0, P 

2: repeat 

3: n+1, P^P 

4- G ^ (P) , U ^ 

5: update {’^,fe,’®'fe,Fc,fe,Ffc,tc,fe,Uc,fe,'Dc.fc,'Ofe}^^^ 

6: P ^ arg^p 

7: ^ minlPcjljLi + iife 

_ 3 _ 

8 : until < en oTt n = nmax 


Proposition 2. Algorithm [T] converges to a stationary point 
of problem A. denoted by (^*, P*, U*, G*), where the cotTe- 
sponding P* is a stationary solution of problem 7?.. 

This can be proved by employing the ideas in CH Theorem 
2]. A sketch of the proof is given as follows. 

Proof of Proposition \2}; The iterations of Algorithm [T] 
monotonically decrease the cost function of A Moreover, the 
feasibility set in constraint (I12cl i is compact, and the mappings 

qMMSE(p) UMMSE(P) 

are continuous. Therefore, the 
iterations converge to a limit point denoted by (P*, U*, G*), 
where P* = arg:4p(U*, G*), U* = and 

G* = G““SE(p*) Furthermore, each of the blocks P*, U* 
and G* satisfies the KKT conditions of its corresponding op¬ 
timization problem, formulated by hxing the other two blocks 
in A. This can be used to show that the point (P*, U*, G*) 
satishes the KKT conditions of problem A. Combining this 
with Proposition [1] completes the proof. ■ 

V. Numerical Results 

We consider a MU-MISO system with Nt = K = 2. 
UncoiTelated channel fading is assumed, where the entries of 
H have a complex Gaussian distribution CM (0, 1). Moreover, 
the noise variance is fixed as cr^ = 1, from which the long¬ 
term SNR is given as SNR = Gaussian CSIT error 
is assumed where the entries of H are generated according 
to the distribution CA/^(0,(Tg). The error variance is given 
as cTg = corresponding to scenarios where the CSIT 

error decays as SNR increases 11, 0. The value of a is 
varied throughout the simulations to represent different CSIT 
accuracies. For each realization H, a channel estimation error 
H is drawn from CM (O, Cg), from which the channel estimate 
is calculated as H = H — H. A channel realization H should 
not be confused with a Monte-Carlo realization While 

the former is unique for a given transmission and not known 



Eig. 1. ASR convergence of Algorithm [T] using 4 different initialization for 
1 randomly generated H, o = 0.6, SNR = 5, 20 and 35 dB, and K = 2. 

to the BS, the latter is part of a sample Hm generated at the 
BS in order to formulate the optimization problem. The size 
of the sample is set to M = 1000 throughout the simulations. 
For a given channel estimate, the mth Monte-Carlo realization 
is obtained as where is drawn from 

the error distribution. 

First, we examine the convergence of Algorithm [T] using 
four different P initializations. For the first initialization (ZF- 
SVD), the precoders of the private messages are initialized as 
ZF-BF constructed using H, while the common precoder is 
constructed by taking the dominant left singular vector of H. 
The second initialization (ZF-e) uses the standard basis vector 
Gi to initialize the common precoder. The third (MF-SVD) 
and fourth (MF-e) initializations apply Matched Beamformers 
(M-BF) instead of the ZF-BF. All initializations use the 
DoF-motivated power allocation from Section III-CI The ASR 
convergence of the proposed algorithm for a = 0.6 and SNRs 
5, 20 and 35 dB, is shown in Figure [T] It is evident that 
the algorithm eventually converges to a limit point regardless 
of the initialization. However, the speed of convergence is 
influenced by the initial state, which also determines the limit 
point, as A is non-convex. The initialization effect becomes 
more visible as SNR grows large. For example, initializing 
the common precoder using SVD enhances the convergence 
at high SNR. In the following results, (MF-SVD) is adopted 
as it provides good overall performance over various channel 
realizations and a wide range of SNRs. 

Next, we consider the Ergodic SR (ESR) performance. 
It is worth noting that the ESR is different to the ASR. 
The latter is the optimization metric defined in Section 
III-Bl which may not necessarily correspond to the actual 
SR achievable at the receivers. The precoders obtained from 
optimizing the ASR yield an achievable SR defined as i? = 
minj -f calculated using the channel 

realization H. Averaging the SR over multiple realizations of 
H yields the ESR defined as Eh {7?}, which is used to capture 
the average performance over multiple channel realizations. In 
the following simulations, the ESR is calculated by averaging 
over 100 channel realizations. The proposed JMB-AWSMSE 




















Fig. 2. MISO-BC and MISO-JMB ESRs. K = 2, and a = 0.6. 



Fig. 3. MISO-BC and MISO-JMB ESRs. K = 2, and a = 0.9. 

scheme is compared to the conventional BC-AWSMSE scheme 
which corresponds to a robust adaptation of the scheme 
proposed in El. Moreover, the base line for conventional 
transmission is taken as naive ZF-BF with Water-Filling (WF), 
i.e. optimization is carried out assuming that the estimate H is 
perfect, and the channel estimation error is not considered. On 
the other hand, we consider a modihed version of the DoF- 
motivated scheme in 13 as a baseline for JMB. In particular, 
the power splitting between the common symbols and the 
private symbols is maintained, while WF is used to allocate the 
power among the private symbols. Furthermore, the common 
precoder is obtained using SVD. The FSRs for a = 0.6, and 
a = 0.9 are shown in Figure |2] and Figure [3 respectively. 
The superiority of all schemes over ZF-BF for the entire 
SNR range is evident. Moreover, JMB-ZF-SVD and JMB- 
AWSMSF achieve the same sum DoF (slope of the curve 
at high SNRs). However, the latter performs better from a 
SR perspective. JMB-AWSMSE and BC-AWSMSE perform 
similarly at low SNRs, where the JMB’s common symbol is 
switched off. The beneht of transmitting a common symbol 
manifests as SNR grows large with a gain exceeding 4 dB 
for a = 0.6, in addition to the DoF gain yielding a faster 
increase-rate. For a = 0.9 which is almost ideal from a DoF 
perspective, the common symbol is not as instrumental as it is 


for lower CSIT qualities. However, ESR and DoF gains can 
still be observed at high SNRs. 

VI. Conclusion 

In this paper, we addressed the problem of ASR max¬ 
imization in a MISO-JMB system with partial CSIT and 
perfect CSIR. The ASR problem was transformed into an 
augmented AWSMSE problem. The AWSMSE problem was 
solved using an AO algorithm which was shown to converge to 
a stationary point of the ASR problem. Numerical simulations 
were employed to demonstrated the benehts of transmitting 
a common symbol in addition to the private symbols. In 
particular, the rate performances of the proposed JMB scheme 
and a state-of-the art linearly precoded MU-MISO scheme 
were compared. At high SNRs, it was shown that the gains 
anticipated by the DoF-based analysis are achieved with an 
enhanced rate performance compared to the DoF-motivated 
design. On the other hand, the proposed scheme converges to 
conventional MU transmission whenever the common symbol 
is not needed, e.g. in the low SNR regime. 
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